Topological Entanglement Entropy with a Twist 
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Defects in topologically ordered models have interesting properties that are reminiscent of the 
anyonic excitations of the models themselves. For example, dislocations in the toric code model 
are known as twists and possess properties that are analogous to Ising anyons. We strengthen this 
analogy by using the topological entanglement entropy as a diagnostic tool to identify properties of 
both defects and excitations in the toric code. Specifically, we show, through explicit calculation, 
that the toric code model including twists and dyon excitations has the same quantum dimensions, 
the same total quantum dimension, and the same fusion rules as an Ising anyon model. 



A fascinating class of many-body quantum systems are 
those that exhibit topological order pQ. Such systems 
are characterized by a gapped ground state manifold, 
with a degeneracy that depends on the boundary condi- 
tions, and have quasi-particle excitations that behave as 
anyons. The degeneracy is robust to local perturbations, 
and therefore such systems are promising candidates for 
storing and manipulating quantum information [2H5]. 

The structure of topologically-ordered systems can be 
further enriched by the use of domain walls or defects, 
across which the quasi-particles can transform nontriv- 
ially. We will be most interested in studying these defects 
in explicit lattice models [Bl-tlOj . and in particular the 
model introduced by Bombin [6]. Essentially the same 
defects have also been studied in Chern-Simons theo- 
ries [TTMl3| . In both of these settings, it has been shown 
that such defects can be viewed as having anyonic-like 
properties that are not associated with the underlying 
model [6l [9l fTTHT3] . More specifically, in a particular 
two-dimensional topologically ordered spin lattice model 
known as the toric code, such a lattice dislocation leads 
to interesting behaviour: the points where dislocations 
terminate, known as twists, interact with the anyons of 
the toric code to reproduce properties, such as fusion 
rules, of nonabelian (specifically, Ising) anyons [6!. This 
is particularly surprising, as all excitations of the toric 
code are abelian anyons. 

In this paper, we further interrogate the analogy be- 
tween twists and Ising anyons by using topological entan- 
glement entropy (TEE) [TJ1 [T5] as a diagnostic. Specifi- 
cally, for a toric code model containing twists, we use the 
TEE to determine: (i) the total quantum dimension of 
the lattice with twists; (ii) the quantum dimensions of the 
objects (quasi-particles and defects) on the lattice, and 
(iii) the quantum dimensions of all of their fusion prod- 
ucts, allowing us to reconstruct the fusion rules for twists 
and excitations. Our results coincide precisely with those 
of the Ising anyon model, lending further support to the 
analogy between the latter and the toric code with twists. 
We note that TEE is a particularly useful quantity in the 
context of the toric code, as both the parent Hamiltonian 



and the modified one with twists are described within the 
stabilizer formalism and so allow the TEE to be calcu- 
lated exactly. 

We first review some relevant properties of topolog- 
ical order and anyons. The low-energy behavior of a 
topologically-ordered system is described by a set of de- 
generate subspaces corresponding to fixed locations of a 
fixed number of elementary particle excitations known 
as anyons. The corresponding Hilbert spaces Ht opo are 
termed "topological" because the different states of the 
anyons cannot be distinguished by local measurements. 
Operations on these spaces can be performed by mov- 
ing anyons along some path and then returning them to 
their initial location, a process known as braiding. The 
topological Hilbert spaces also have a map corresponding 
to bringing two anyons together to create a new one, a 
process known as fusion. The dimension of 'Htopo for a 
system with N type-a anyons in the limit of large N is 
given by dim 'Htopo ~ {d a ) /T> 2 . Here, d a is the quantum 
dimension of a type-a anyon, and T> = y^J^x wnere 
the summation is over all the anyon types of the given 
model, is the total quantum dimension. Interestingly, d a 
need not be an integer, which underlines the fact that 
"Htopo need not have a tensor product decomposition into 
subspaces associated with individual anyons, but rather, 
the topological states are a collective property of all the 
anyons. 

Two anyon models are relevant to our study: Z 2 
abelian anyons, and Ising anyons (see, e.g., Ref. [To]h 
The elementary excitations of a Z 2 abelian anyon model 
are the electric charge, e, the magnetic vortex, m, the 
dyon, e and the vacuum, 1. These anyons have quantum 
dimension d a = 1 for all anyons a, and therefore 2?z 2 = 2. 
A model possessing Ising anyons includes as elementary 
excitations: an Ising anyon er, a fermion ip, and the vac- 
uum 1, where d\ = = 1 and d a = a/2. These anyons 
are nonabelian, as when two a anyons fuse, there are two 
possible fusion outcomes 1 and ip. Like the Z 2 quantum 
double, this model has T>i smg — 2. 

Entropy is a powerful tool for probing topological 
properties of a model. For the ground states of quan- 
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FIG. 1: (a) Re gion s used to calculate topological entangle- 
ment entropy. |(b)| Regions used to probe the quantum di- 
mension of an anyon on region D\. 



turn many-body systems with local Hamiltonians the von 
Neumann entropy Sa of the reduced state on a region A 
is expected to satisfy an area law [IT] [T^ . Restricting to 
two-dimensional systems, the von Neumann entropy of a 
region A of a topologically ordered system is expected to 
have a universal correction 7 to Sa such that [121 US] 



aL A - n A j, 



(1) 



where La is the length of the boundary of region A 
(which should be smooth), ua is the number of discon- 
nected boundaries enclosing region A and a is a constant 
that depends on the microscopic properties of the Hamil- 
tonian. 

For a given topological model, the total quantum di- 
mension is related to 7 via 7 = log V [TH [T5] . As such, 
the von Neumann entropy can be used to probe the any- 
onic properties of the system. It is useful to eliminate 
the the boundary contribution of Eq. ([I]) by computing a 
linear combination of the von Neumann entropies of over- 
lapping regions as in Fig. 1(a) to obtain the topological 
entanglement entropy 

Stopo = Sa + Sb + Sc — Sab — Sbc — Sca + Sabc , (2) 

where Sr is the von Neumann entropy of region R [141 
[T5] . The TEE is a function solely of the total quantum 
dimension of the model, S\opo = — 7 = — log P. 

Additional information about the anyonic properties of 
the model can be revealed by defining analogous entropic 
quantities using regions with noncontractable bound- 



aries [in], as in Fig. 1(b) By using excited states rather 
than the ground state, our entropic calculations will de- 
pend explicitly on the number of anyons present and their 
placement within the regions. Consider a state describ- 
ing an anyon of type a on a simply-connected region, 
labeled A; the presence of the anyon gives a correction to 
the von Neumann entropy as 5.4(a) = aLA — log(f/rf a ) 
[T4l [20] . Again, the boundary contribution can be re- 
moved with a suitable combination of the entropies on 
a nontrivial partition; a suitable choice is an annular 
region shown in Fig. l(b)| defining a disconnected re- 
gion D = D\ U Z?2- For a state describing a pair of 




FIG. 2: Pauli X, Y and Z operators are shown in blue, red 
and green circles respectively, (a) A four-body interaction 
term for Kitaev's toric code Hamiltonian. (b) and (c) Strings 
of Pauli operators which respectively generate electric and 
magnetic charged pairs at the end points. Electric(magnetic) 
strings can only be generated on grey(white) plaquettes. (d) 
The five-body interaction term for the twist plaquette. (e) A 
Pauli string operator that cross a domain wall, (f) A logical 
operator in the form of a dyon chain terminating at two twists. 



anyons of type a that fuse to the vacuum, with one 
in region D\ and the other in region D2, the entropy 

S'aim(a) = Sa + Sb + Sc — Sab — Sbc — Sca + Sabc 
gives a value 



S aDn (a) = -2logV/d a , 



(3) 



where the additional factor of 2 can be seen as arising 
from the boundary connectivity as in Eq. Q. 

In this work we propose to use Sarin to infer fusion 
rules for both anyons and twists. For this purpose we 
must determine how S ann behaves in the situation where 
there are multiple anyons in region D\. If 'Htopo is i n 
a state where the anyons ai, 02, • ■ ■ in D\ fuse to a defi- 
nite fusion product X\, then S aml (ai, (Z2, ■ ■ ■) = S ann (xi). 
However, in general, anyons can fuse through multiple 
channels, with the probability (amplitude) of a given fu- 
sion outcome determined by the state encoded in Ht opo - 
Considering only states where the connected region ABC 
as well as the entire lattice fuse to the vacuum, the out- 
come X\ determines explicitly X2, the fusion product of 
the anyons on Di. In this important special case, we can 
use results from Ref. [20] to compute S'ann and we find 
the explicit formula 

S an „({x}) = -21og2? -Y, p ^°dPJ{d Xl d X2 )\ , (4) 



where the sum is over all fusion outcomes x = (xx,X2), 
each of which occur with probability P x . This equation 
provides the principal diagonostic we will use in this pa- 
per to characterise the anyonic properties of a model's 
quasi-particles and defects. As we will show, the TEE 
of a model with twists using annular regions with a non- 
contractable boundary will yield identical results to the 
analogous entropic calculation using excited states of a 
model with Ising anyons in place of twists. 
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The specific model we consider here is the Z2 quan- 
tum double, also known as Kitaev's toric code model [3], 
and we focus on a formulation due to Wen [2TJ [22] . 
On a square lattice with spins on the vertices, labeled 
(j,k), the model is described by the Hamiltonian H = 

~ J2j.k w j,k with Wj tk = Xj t kZj + i t kZj t k+iXj + i t k+i- We 
note that H is a stabilizer Hamiltonian, and all term com- 
mute. The low energy excitations of this model are Z2 
anyons. Non-trivial domain walls can be introduced in 
this model by shifting the geometry of the lattice along a 
line, such that we have even(odd) plaquettes adjacent to 
other even(odd) plaquettes IB] as shown in Fig. [2] The 
interaction terms of the Hamiltonian along the disloca- 
tion then change to correspond with the new geometry, 
but remain locally equivalent to the old interaction terms 
except at the ends of the dislocation. The five body in- 
teraction terms where dislocations terminate, such as in 
Fig. [2j are known as twists [6] . We note that the insertion 
of twists into the toric code model changes the ground 
state degeneracy, as can be seen from the fact that the 
dyon chain stabilizers shown in Fig. [2^f ) commute with 
the Hamiltonian. 

The presence of twists in the toric code model pre- 
vents us from globally defining the electric or magnetic 
change of an anyon: an e(m) anyon that crosses a do- 
main wall will transform into an m(e) anyon, as shown 
in Fig. [2je) , although the dyons remain well defined. In 
addition, the twists at the end points of these disloca- 
tions can be used to reproduce some of the fusion rules 
of the Ising anyon model [6]. This analogy is achieved 
by considering only the vacuum, toric code dyons e, and 
one type (handedness) of twist. By treating e dyons as 
the fermions (V'-type anyons) of an Ising anyon model, 
and the twists as a-type anyons, one can reproduce the 
fusion properties of the Ising anyon model. In this sense, 
twists exhibit anyon-like behaviour. Here we use TEE as 
a diagnostic tool to further interrogate this analogy. 

We can exploit the stabilizer formalism to calculate von 
Neumann entropies for regions of the vacuum and excited 
states of these topological models, using the method of 
Ref. [23]- This method makes use of a canonical form 
of a generating set of the stabilizer group defined by a 
bipartition of the lattice into regions A and A. The sta- 
bilizer generators can be grouped into those entirely in 
A, those entirely in A and those that have support on 
both partitions. For the set of stabilizer generators that 
cross the bipartition, a canonical form can be constructed 
such that the local support of the generators Ki on A, 
denoted by Ki, either commutes with all other Kj, or 
anti-commutes with only one other such element -KV. 
These pairs of stabilizer generators that anti-commute 
when restricted to A describe a pair of maximally entan- 
gled qubits sitting across the boundary. Counting these 
pairs gives the von Neumann entropy Sa- Applying this 
method to the toric code vacuum [24) gives Sa — hLA—1, 
and therefore Stopo = ~ 1 and T> — 2 as expected for a 




FIG. 3: The boundary dA shown in red of a simply connected 
region enclosing a twist on a lattice. We number the cut 
plaquettes with a single index shown in blue, here, La = 13. 
We also show a logical operator deformed outside region A. 



model with Z2 anyons. (Here we choose the length of the 
boundary La to correspond to the number of plaquettes 
that cross the boundary.) We note that the following 
results are independent of the topology of the lattice, 
because homologically non-trivial stabilizers of a topo- 
logically non-trivial model can always be deformed away 
from the regions considered in the present entropy cal- 
culations [24]. With this in mind, we can obtain the 
same results for the following calculations regardless of 
the boundary conditions of the lattice. 

We first consider a vacuum state of a toric code model 
with some number of twists, and calculate the von Neu- 
mann entropy for a simply connected region enclosing 
one twist as shown in Fig. [3] We will show that Sa — 
\La — §, as expected for a region of an Ising anyon model 
containing an cr-type Ising anyon with d a = To 
do this calculation, we must identify a generating set of 
stabilizers for a vacuum state, in the canonical form de- 
scribed above. First, although there exist stabilizers that 
take the form of dyon strings with ends terminating at 
twists, as shown in Fig. [2]jf), we can multiply by local 
stabilizers to deform such strings such that they have no 
support on A, as shown in Fig. [3j For plaquette opera- 
tors Wi that cross the partition dA, we define Wi to be 
the support of these operators on A, labelling them with 
a single index 1 < I < La as in Fig. [3] The Wi opera- 
tors are not in canonical form, as each Wi anti-commutes 
with both Wi-i and To achieve canonical form, we 

replace one generator, W m , with the generator Y\ gA W, 
where the product is over all stabilizers with support on 
the boundary of A; the support of this operator on A 
commutes with all Wj. (Note that, in the untwisted toric 
code, one can construct two such operators using the nat- 
ural partition of toric code stabilizers into two commuting 
subsets. With the presence of a twist, however, there is 
only a single such operator.) The remaining Wi can be 
transformed such that they all pairwise anti-commute, 
giving the desired canonical form. Counting the pairs of 
locally anti-commuting generators then gives the stated 
result Sa — \La — \ - 
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FIG. 4: |(a)| The annular topological entanglement entropy 
regions with one twist on region D\ and one twist on region 
Di- (b) The lattice configuration where two twists lie on Di 
and two lie on D2, in an eigenstate of Z (defined as a dyon 

A configuration with four 



chain stabilizer as in Fig. j2^f)). 
twists in an eigenstate of X. 



This method allows us to calculate the effect of twists 
on the TEE of the vacuum, using Eq. ^ and regions de- 
fined as in Fig. |l(a)| We find that for a vacuum state of a 
model with any number of twists, the TEE is Stopo = — 1 
independent of the locations of the twists. This outcome 
corresponds to a model with total quantum dimension 
D = 2, consistent with models with either Ising or 7L 2 
abelian anyons. 

We next show that, using the annular regions with non- 
connected boundaries as in Fig. [T(b)1 which can probe the 
properties of anyons, the entropy Sam^l) for annular re- 
gion of the vacuum state of a model containing twists is 
equivalent to the entropy S ann (a) of a state of an Ising 
anyon model with cr-type Ising anyons. We use a parti- 



tioning as in Fig. 1(b) with a twist on region D\ and an- 
other on D 2 , as shown in Fig. 4(a) We achieve canonical 
form for region ABC using the method described above, 
but now with disconnected boundaries we must consider 
the canonical form of the stabilizer generators for each 
boundary separately. As above, any dyon string stabi- 
lizer that crosses region ABC can be deformed entirely 
onto region D 2 . Calculating the entropy for the various 
regions, we find Sabc — \Labc — E with all other re- 
gions giving S R = \L R - 1, and so S top o = -1- We 
compare this result to a calculation of Sarin (c) for a state 
of an Ising anyon model with a cr-type Ising anyon in each 
region D\ and D 2 , using Eq. ^ , and obtain an identical 
result, 5 , ann (cr) = —1. Thus, from an entropic perspec- 
tive, the twists appear as cr-type Ising anyons. A similar 
calculation involving a e quasi-particle on region D\ and 
another on region D 2 gives S'ann(e) = <S'ann('/') — 1 as ex- 
pected from a Ising fermion ip in the Ising anyon model. 

Finally, we use TEE to probe the fusion rules associ- 
ated with quasi-particles and twists. We show that pairs 
of twists fuse according to the more sophisticated non- 
abelian fusion rules of two cr-type Ising anyons, and mod- 
ify the TEE according to our Eq. @ in a way that de- 
pends on the specific fusion channel. We place two twists 
on region Di and two twists on region D 2 . The channel 
by which the two twists on region D\ fuse together is 
dependent on the logical state, as characterised by the 



logical operators that are given by dyon strings between 
twists. 

We first consider the lattice prepared in a n eige nstate 
of the logical-if operator Z shown in Fig. 4(b) The 
+1 eigenstate of Z corresponds to two cr-type anyons 
on region D\ fusing to the vacuum 1 with certainty, 
whereas the —1 eigenstate fuses to the fermion ip with 
certainty. In either case, we find Sa, nn ({(a, u) — > 1}) = 
•Sarin ({ (c, o) -> ip}) — —2. Comparing this with Eq. Q, 
we have the result we anticipate from the Ising anyon 
model. 

We next consider the case where the lattice is prepared 
in an eigenstate of the logical- X operator X shown in 
This is analogous to the Ising anyon model 



Fig. 4(c) 



where one cr-type anyon on region Di and one cr-type 
anyon on region D 2 will fuse to give one fusion outcome 
with certainty. However, fusing the two particles on re- 
gion D\ can produce the vacuum 1 or a fermion ip with 



equal probability P(i.i) = P{ 



The details of 



the TEE calculation arc more involved in this case, as 
the operator X cannot be deformed away from region 
ABC. We therefore require a different method to find a 
canonical form for region ABC. Denoting the two sim- 
ply connected segments of the boundary by dABCp, for 
P = 1,2, we rewrite the indices of the plaquettes cut 
these boundaries as W^;, with 1 < I < Lp number- 
ing the generators around the boundary. We deform X 
such that it only has common support with the genera- 
tors Wi i l 1 -i, Wi t Lx, W 2 ,L 2 -it and W 2 x 2 on t ne bound- 
aries. To achieve canonical form, we replace the genera- 
tors Wl, Ll -l -> UdABCx W l,h W 2 ,L 2 -1 -> UdABC 2 W 2,h 
and W hLl UdABCuk^venWhkUdABC^ven^U 

such that their local support commutes with the lo- 
cal support all the other generators. The remaining 
cut generators, including X, are then transformed such 
that their local supports pairwise anti-commute, giv- 
ing the required canonical form. This gives the result 
Sabc = \Labc — 1- The other regions used fuse to vac- 
uum, giving S ann ({(a,a) — > 1 or ip}) = —1, which from 
Eq. Q is the result we expect from the analogous Ising 
anyon case. 

To conclude, we have shown that we obtain the same 
topological entanglement entropy results we expect for 
the Ising anyon model using the toric code lattice with 
twists. Specifically, we have shown that the lattice ad- 
mits the same total quantum dimension, and the objects 
on the lattice and their fusion products give the same 
quantum dimensions as the Ising anyon model, further 
strengthening the analogy highlighted in [B]. 

To complete the analogy, one might seek to identify the 
braiding properties of twists, by using the techniques of 
Ref. |19j for example, and compare them to Ising anyons. 
However, as noted by Bombin j6j, braiding of twists are 
not expected to be well defined and we have not consid- 
ered them in this work. We note that it may be possi- 
ble to define robust braiding operations up to an overall 
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phase [TU1 H2J US] , and in particular one might seek to 
use techniques from Ref. [13] to check explicitly that the 
twists in this model have the expected braiding proper- 
ties. We leave this for future work. 
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